Polynomiality of unpolarized off-forward distribution functions 
and tfie D-teim in tfie cfiiral quark-soliton model 

p. Schweitzer", S. Boffi"'^ M. RadicP''' 

OA ' " Dipartimento di Fisica Nucleare e Teorica, Universitd degli Studi di Pavia, 1-27100 Pavia, Italy 



o 
o 

(N 
> 

o 
m 
(N 



Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, 1-27100 Pavia, Italy 

September 30, 2002 
Abstract 



Mellin moments of ofF-forward distribution functions are even polynomials of the skewedness 
parameter ^. This constraint, called polynomiality property, follows from Lorcntz- and time- 
^ ' reversal invariance. We prove that the unpolarized off-forward distribution functions in the 

, chiral quark-soliton model satisfy the polynomiality property. The proof is an important con- 

■ tribution to the demonstration that the description of off-forward distribution functions in the 

(N ■ model is consistent. As a byproduct of the proof we derive explicit model expressions for mo- 

' ments of the D-tevm and compute the first coefficient in the Gegenbauer expansion for this term. 

O ■ PACS: 13.60.Hb, 12.38.Lg, 12.39.Ki, 14.20.Dh 

1 Introduction 

Oh! 

^ I Off-forward distribution functions (OFDFs) of partons in the nucleon - or equivalently double, 
skewed or non-diagonal parton distribution functions [1, 2, 3, 4] - are a promising source of new 
information on the internal nucleon structure (see Refs. [5, 6, 7] for recent reviews). OFDFs 
^ I enter the description of hard exclusive processes, such as hard meson production [4] and deeply 
' virtual Compton scattering [8]. Only recently first data on the latter process became available 
[9, 10, 11, 12], and so our present understanding of non-perturbative aspects of OFDFs has to rely 
on models. 

After the first estimates in the framework of the MIT bag model [13], studies in the chiral 
quark-soliton model (xQSM) of the nucleon have been presented [14, 15]. The xQSM has been 
derived from the instanton model of the QCD vacuum and provides a relativistic quantum field- 
theoretical description of the nucleon. Without any adjustable parameters it describes a large 
variety of nucleonic properties, like form factors [16, 17] and (anti)quark distribution functions 
[18, 19, 20] typically within (10 - 30)%. 

However, the reliability of a model is not only based on its phenomenological success. From 
a theoretical point of view, it is important to demonstrate the internal consistency of the model. 
One of the virtues of the xQSM, which is due to its field-theoretical character, is the possibility to 
provide analytical proofs which explicitly demonstrate the consistency of the model. In Ref. [18] it 
has been proven that the quark and antiquark distribution functions computed in the model satisfy 
all general requirements, such as sum rules, positivity and inequalities. With the same rigour it 
has been shown in Refs. [14, 15] that the xQSM expressions for OFDFs reduce to usual parton 
distributions in the forward limit, and that their first moments yield form factors. 

In this paper we will present a further check of consistency for the xQSM, and give for the 
unpolarized OFDFs an explicit proof of polynomiality, i.e. of the property that the m*^ Mellin 
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moment of an OFDF is an even polynomial in the skewedness variable of degree less than or 
equal to m at the (unphysical) value t = of the Mandelstam variable ^. The polynomiality 
property follows from Lorentz- and time-reversal invariance of strong interactions. The proof of 
polynomiality is the main result of this paper. 

An elegant way to satisfy the polynomiality condition is to use the double distributions (DDs) of 
Refs. [1, 2]. However, if assumed to be regular functions, DDs yield a vanishing highest power of ^ 
for singlet moments of unpolarized OFDFs. This means that DDs yield an incomplete description, 
which can be completed by introducing the so-called D-tevm [21]. Moments of the D-tevm are 
uniquely defined in terms of the coefficients of the highest power of ^ [21, 22, 23]. In Ref. [24] the 
D-term has been calculated numerically from moments of the OFDFs evaluated in the model for 
physical values of t. The results obtained here enable us to evaluate in the model, in principle, all 
Mellin moments of the D-term directly at the unphysical point t = 0. We will demonstrate this for 
the lowest non- vanishing moment. 

The paper is organized as follows. In Sec. 2 a brief introduction to the xQSM is given. In Sec. 3 
the model expressions for the unpolarized OFDFs are presented. The proof of polynomiality is 
given in Sees. 4 and 5 for the non-spin-fiip and the spin-flip OFDF. Sec. 6 is devoted to the D-term. 
Finally, in Sec. 7 we summarize our results and conclude. Some technical details required for the 
proofs are given in the Appendices. 



2 The chiral quark-soliton model (xQSM) 

The xQSM [25] is essentially based on the principles of chiral symmetry breaking and the limit of 
a large number of colours Nc- The effective chiral relativistic quantum field theory, which underlies 
the xQSM, is formulated in terms of quark, antiquark [jp, 'ijj) and Goldstone boson (pion tt", 
a = 1, 2, 3) degrees of freedom. It is given by the partition function [26, 27, 28] 

Zeff = JvipVij; VU exp J d'^x ip {i^ - M )V'^ , U^^ = e^^sr^Tr" 

In Eq. (1) M is the dynamical quark mass due to spontaneous breakdown of chiral symmetry, which 
is in general momentum dependent, and U = exp(ir"7r") denotes the SU{2) chiral pion field. The 
effective theory in Eq. (1) contains the Wess-Zumino term and the four-derivative Gasser-Leutwyler 
terms with correct coefficients. It has been derived from the instanton model of the QCD vacuum 
[28, 29], which provides a mechanism of chiral symmetry breaking. The effective theory (1) is valid 
at low energies below 600 MeV, a scale set by the inverse of the average instanton size. 

One application of the effective theory (1) is the chiral quark-soliton model (xQSM) of baryons 
[25]. In large- TVc limit the path integral over pion field configurations in Eq. (1) can be solved in 
the saddle-point approximation. The large- A^c limit is known to be a good theoretical guideline. In 
this limit the nucleon can be viewed as a classical soliton of the pion field [30]. The xQSM provides 
a realization of this idea. In the leading order of the large- A^c limit the pion field is static, and one 
can determine the spectrum of the effective one-particle Hamiltonian of the theory (1) 

i?eff|n) = En\n) , H,s = -il^j'^dk + -/^MU''' . (2) 

The spectrum consists of an upper and a lower Dirac continuum, which are distorted by the pion 
field as compared to continua of the free Dirac-Hamiltonian 

//()|/'o) = E,,,\no) , ^0 = -n°7*^5fe + 7°M , (3) 
^Here we use the notation of Ref. [3]. 
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and of a discrete bound state level of energy E\cy, for a strong enough pion field of unity winding 
number. By occupying the discrete level and the states of lower continuum each by Nc quarks in 
an anti-symmetric colour state, one obtains a state with unity baryon number. The soliton energy 
Egoi is a functional of the pion field, 



Minimization of -Egoji/] determines the self consistent solitonic pion field Uc- This procedure is 
performed for symmetry reasons in the so-called hedgehog ansatz 7r"(x) = ^ ^'(Ixj) with the radial 
(soliton profile) function P(|x|). The nucleon mass Mn is given by £'soi[C4]- Quantum numbers of 
the baryon like momentum, spin and isospin are described by considering zero modes of the soliton. 
Corrections in 1/Nc can be included by considering time dependent pion field configurations. The 
results of the xQSM respect all general counting rules of the large- A^c phenomenology. 

For the following it is important to note that the effective Hamiltonian Hf.g commutes with the 
parity operator n and the grand-spin operator K, defined as the sum of quark angular momentum 
and isospin. Thus ^eff ) ft, and form a maximal set of commuting operators, and the single 
quark states |n) in Eq. (2) are characterized by the quantum numbers parity tt, K and M 



The xQSM allows to evaluate in a parameter-free way nucleon matrix elements of QCD quark 
bilinear operators as 



where for zi ^ Z2 the insertion of the gauge link is understood on the LHS. The dots in Eq. (6) 
denote terms subleading in the 1/Nc expansion, which we will not need here. In Eq. (6) F is some 
Dirac and flavour matrix, A^^^y, a constant depending on F and the spin and flavour quantum 
numbers of the nucleon state \N) = {83, T^), and $n(x) = (x|n) are the coordinate-space wave- 
functions of the single quark states \n) defined in Eqs. (2, 5). The sum in Eq. (6) goes over occupied 
levels n (i.e. n with En < Eiev), and vacuum subtraction is implied for En < E\ev as in Eq. (4). 

Many static nucleonic observables - like magnetic moments, electric polarizabilities, axial prop- 
erties, etc. - have been computed in the xQSM in the way sketched in Eq. (6). The results were 
found in good agreement with data (see Ref. [17] for a review). In particular the model describes 
data on electromagnetic form factors up to \t\ ~ 1 GeV^ within (10-30)% [16, 17]. In Ref. [18] it has 
been demonstrated that the model can be applied to the description of twist-2 quark and anti-quark 
distribution functions of the nucleon. The consistency of the approach has been shown by giving 
proofs that the model expressions satisfy all general requirements of QCD [18]. The distribution 
functions computed in the xQSM [18, 19, 20] refer to a low normalization scale of around 600 MeV, 
and agree with parameterizations performed at comparably low scales [31] within (10-30)%. 

The success of the xQSM in the parameter-free description of - among others - form factors and 
distribution functions encourages confidence into the predictions for OFDFs made in Refs. [14, 15]. 

3 Off-forward distribution functions in the xQSM 

The unpolarized quark off-forward distribution functions are defined as [3] 




(4) 



n) = \En, TT, K, M) . 



(5) 




n,occ 




= H\x, t) U{V', S'^) fi C/(P, S3) + Eiix, ^, t) f7(P', S'^) 




U{P,Ss) + ... (7) 
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where for brevity the gauge hnk is omitted and the scale dependence not indicated. The dots 
denote higher-twist contributions. The hght-hke vector n^, the four-momentum transfer A'^, the 
skewedness parameter ^ and the Mandelstam variable t are defined as 

= , n{P' + P)=2, A" = {P' - , nA = -2^ , t = A'^ . (8) 

In Eq. (7) X e [—1, 1] with the understanding that for negative x Eq. (7) describes minus the OFDF 
of the antiquark. When evaluating the above expressions in the xQSM in the large-A^c limit, one 
has to note that Mn = C(-^c) is much larger than spatial components of the momentum transfer 
I A* I = 0{N^) which in turn is much larger than the zero-component of the momentum transfer 
|A°| = 0{N~^). The variables x and ^ are of 0{N~^). Choosing the 3-axis for the light-cone space 
direction, we have in the "large- AT^, kinematics" 

n^ = {l,0,0,-l)/M^ , t = -A\ e = -AV(2MN). (9) 

Different flavour combinations of the OFDFs exhibit different behaviour in the large-iVc limit [7] 

(iJ" + H''){x, e, t) = N^fiNcX, iVeC, t) , {E^ - E''){x, ^, t) = N^fiN^x, TV,^, t) , (10) 
{H^ - H^){x,^,t) = Ncf{N,x,NcC,t) , {E^ + E'^){x,C,t) = N^f{N,x,N,^,t) . (11) 

The functions f{u, v, t) in Eqs. (10, 11) are stable in the large- A'^c limit for fixed values of the 0{N^) 
variables u = NcX, v = and t, and of course different for the different OFDFs. 

The model expressions for the leading OFDFs, Eq. (10), have been derived in Ref. [14] and read 

{H^ + H'')ix,^,t) = M^A^e/d^Xe^^^ ^ j ^ ^iz%xM^-E„) 

X $;(X + ^e^) (1 + 7V) ^>n(X - ^e^) , (12) 

{E--E'){x,i,t) = ^^/d^Xe^^^^ /^,.A^M.-^n) 

X $;(X + ^e^) (1 + 7V)(T X A)3$„(X - ^e^) . (13) 

Here is the unit vector of the 3-axis, which is singled out by the choices made in Eq. (9), A"*" 
denotes the part of A transverse to 3-axis, and (r x A)^ = gS'^^V^A'. 

In Ref. [14] it has been demonstrated that in the forward limit {H^ + H'^){x,^,t), Eq. (12), 
reduces to the model expression for unpolarized isoscalar distribution function 

lhn(if- + H''){x,^,t) = if^ + ft){x) (14) 
t->o 

and that the model expressions (12, 13) are correctly normalized to the corresponding electromag- 
netic form factors 

1 1 
jdx {H^ + H''){x,^,t) = (F« + Ffm , Jdx {E^ - E''){x,^,t) = {F^ - Fi){t) . (15) 
-1 -1 

The next and maybe most stringent check of the model expressions is the demonstration that 
Eqs. (12, 13) fulfil the polynomiality condition. As a consequence of Lorentz- and time-reversal 
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invaxiance the rn"^ Mcllin moment M^"*) t) of an OFDF is an even polynomial in ^ of degree less 
than or equal to m at i = 

dx if-^lx, ^, 0) = hi + /.I ('"^C^ + . . . + (16) 
^ I s for m odd, 

^ L ^m-i s odd. 

Due to the spin ^ of the nucleon the highest coefficients of even moments of H'^{x, ^, t) and -E*(x, ^, t) 
are related to each other by 

h'Jr^ = -et'^ . (18) 

The point t = is unphysical and can be reached only by means of analytical continuation. In 
Ref. [14] (i?" + H'^){x,i,t) and (E" - E'^){x,i,t) have been computed as functions of (physical 
values of) x, ^ and t, and the polynomiality conditions, Eqs. (16, 17), have been checked by taking 
(numerically) moments M^^(^,t), M^^\(,,t) and extrapolating (numerically) to the unphysical 
point t = [24]. It is clear in this way one cannot prove strictly polynomiality. This will be done 
in the next two sections. 

4 Proof of polynomiality for {H^ + H^)(x, ^,t) 

The m^^ moment of (i?" + H'^){x,^,t) in Eq. (12) - which we will refer to as M^^(^,t) - reads 
(see App. A) 

MP{i,t) = E {"^k ') ^ E (J) ("1(1 + 7V) (p^)^ exp(zAX) ip^r^ln) , 

(19) 

where p^, X mean the 3^'^ component of the free-momentum operator and the position operator, 
respectively. The next step is to take the limit t ^ with ^ 7^ in Mf^^''\^,t). Eq. (19) shows 
that M^\^,t) = M^\a.), i.e. A determines entirely the dependence of the moments on ^ and 
t according to Eq. (9). Moreover, A appears only in the operator exp(zAX). The order of the 
operations of taking the limit t ^ and taking matrix elements in Eq. (19) can be interchanged. 
Taking the limit t — (understood as analytical continuation) on exp(zAX) we obtain (see App. B) 

hm exp(iAX) = V (-^^^^nIX] )'° Pi,(cos^) . (20) 

analytical fcg! 
continuation le—0 

Inserting Eq. (20) into Eq. (19) and simplifying the result by using symmetries of the model (see 
App. C) we obtain 

Mr«.o) ^ ^ E zn')^' t 

n.occ fe=o ^ ie=0 

even 

k 

xE(j) {n\{lW)HP')m'' PlA^ose){p^f-^\n) , (21) 

j=0 
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where (7^7^)'^, equal to ^y^^^ for odd and unity for even k, is introduced to simplify the notation. 
Note that in the expression for M^^\$^,Q) only even powers of ^ appear. Odd powers of vanish 
due to parity transformation (see App. C). What remains to be shown is that this series in is a 
polynomial of degree less than or equal to m. 

The operators in the matrix elements in Eq. (21) transform as irreducible tensor operators 
of rank L and M = under simultaneous rotations in space and isospin-space. More precisely |X| 
(and any power of it) is rank zero, 7'^7^ and fr' are rank 1, P;^(cos^) is rank 1^ (see App. D). The 
product [^]m2 irreducible tensor operators and [Fi^m^ can be decomposed into 

a sum of (certain new) irreducible tensor operators according to 

m]M\[^2]fe= ^'m with f^^C(LM;LiMi,L2M2)[fi]^jf2]^^, (22) 

L=\Li~L2\ 

where the C(LM; L\M\, L2M2) denote Clebsch-Gordan coefficients. The product of three or more 
irreducible tensor operators can be decomposed in an analogous way. 

For even k in Eq. (21), there is a product of A; + 1 irreducible tensor operators, namely the 
rank-1 operator , which appears fc-tinics, and the rank-^g operator 1X1'"= P;_.(cos ^). For odd k 
there is the further operator 7*^7'^, i.e. there are altogether k-\-2 irreducible tensor operators. Upon 
successive application of Eq. (22) these A; + 1 or A; + 2 irreducible tensor operators yield a series of 
certain new irreducible tensor operators [2]?]^ with even ranks L ranging from zero to 

_ ( k + le + l if A; odd , . 

'"^ " U + if even. ^^"^^ 

Note that only even ranks appear since all involved operators have the "magnetic quantum number" 
M = 0. The operators in Eq. (21) can thus be decomposed as 

^max 

(7%^)^ ip'y\±\'^ PiA^^osO) {f f-^ =J2^L [Th]^ ■ (24) 

L=0 

Writing out explicitly the dependence of the single quark states |n) on all quantum numbers, 
Eq. (5), we have 

E ("I (tV)' ip'ym'' mcose) if)'-'\n) 

n,occ 

-^max 

= E E E «L i^n, 7T, K, M\ [fH]li\En, TT, K, M) . (25) 

n,occ K.M L = 

On the RHS of Eq. (25) we recognize traces (sums over matrix elements diagonal in K, M) of 
irreducible tensor operators, and the latter vanish unless the operator has rank zero [32]. As a 
consequence, the sum over L in Eq. (25) receives a non-zero contribution from \Th^ only. 

Consider the case when k is odd. The situation for L = in which If, takes its largest possible 
value Vf^^^ happens when 7''7'^ and the (in whatever way ordered) k p^-operators combine to an 
operator of rank A; + 1. Then Z™^ = A; + 1. If A; is even, (7^7^ )'= = 1 does not contribute and then 
^max _ ^ g^^ £qj. ^ given k, the infinite series in Zg in Eq. (21) actually stops at an Z™^ which is a 
function of k 

;max/T,^ _ \ k + l for A; odd /^^.x 
^^>-\k for it even. 

Inserting the result (26) into Eq. (21) yields 

even 
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k 

X E (^1(7%')' |X|'^ PUcose) (p'^f-^n) . (27) 
j=o 

We observe that Mj7^(^,0) is a polynomial in even powers of ^, the highest power being given by 

IT-'ik) =1"^ , form even 

m — 1 for m odd, 

which completes the proof of polynomiality, Eq. (16), for {H^ + H'^){x,^,t). 

5 Proof of polynomiality for (E"^ — E^){x,^,t) 

The model expression for the m*^ moment of {E'^ — E'^){x, ^, t) in Eq. (13) is given by (see App. A) 

3Mn „,occ fe=0 i=0 

X (n|(l + tV) (P^)^(T X A)3 ^^1^^ (ff-^\n) . (29) 

To perform the limit t — in Mg"^(^,t) we have to continue analytically zA"*- exp(iAX)/(A-'-)^ 
to t = 0. The result reads (see App. B) 



zA^exp(iAX) ^ (-2z^Mn)'«-2 
lim j = > — 

analytical (A )2 , te ' 



continuation 



(30) 



Inserting Eq. (30) into Eq. (29) and simplifying the result by means of symmetries (see App. C) 
we obtain 

^^Wn „,occ jfc=o ^ ^ fe- 



/e even 

k 

" ■ '(rxp)^ |X|'«P,,(cos0) 



xE(J)H(7%Y+^(p^)^'^ 



iff-^. (31) 



Note that the contribution of Ze = 1 in Eq. (30) is singular at ^ = 0, but it does not contribute in 
Eq. (31), such that M^"^\^,0) is well defined for all values of ^. We observe that M^™-'(^, 0) is a 
series in even powers of ^. Thus, what remains to be done again, is to demonstrate that this series 
is actually a polynomial of degree less than or equal to m according to Eq. (17). 

For this purpose we note that also (r x p)^ is an irreducible tensor operator of rank 1 with 
M = (see App. D), such that in Eq. (31) we deal with matrix elements of a product of irreducible 
tensor operators. As argued in the previous section, these operators can be decomposed into a 
sum of certain new irreducible tensor operators [Te]m with M = and (even) ranks L where 
< L < Lmax with 

_r/c + /e + l if/c odd . s 

'"^"U + Ze + 2 if even. ^'^^> 

The operator in Eq. (31) contains a commutator of (r x p)^ and |X|''=P;^ (cos ^). Focusing first on 
one ordering of these operators the decomposition into a sum of irreducible tensor operators, which 
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we will call [T^^^Jm, looks as follows 



L=0 

where the coefficients b^^^ and the operators [Tg^^]o follow from the successive application of the 
tensor operator multiplication rule in Eq. (22). However, only the [TfeJo contributes to the sum 
over single quark states in Eq. (31). 

There are many ways to construct a rank zero operator out of ^ {p^y , (t x p)^, 

|X|'= Pi^{cos6) and {p^)^^K The relevant one (for the proof) is the one which allows 1^ to take its 
maximal value This case happens when the ranks of (7^7^)'^'''^, {p^y , (t x p)'^ and {p^)^~^ 

add up to A; + 2 for even A;, or to A; + 1 for odd k. Then lf^{k) takes the value 

;max/. N _ I k + l for A; odd 
^''^"lA: + 2 for A; even. 

Considering the other ordering of (r x p)^ and |X|'«^P;^(cos ^) in the commutator in Eq. (31), one 
of course arrives at different coefficients b^l^ and operators [r^^]^ in Eq. (33), but at the same 
conclusion of Eq. (34). Applying the result (34) to Eq. (31) we obtain 

2JV, „ Im-WK'-'-'' 'Vi" (-2'{-U.)'--' 



^Jyl^ n,occ jfc=o ' z te- 



xE(j)H(7V)'+'(p'y^ (Txp)3, |X|'«P,^(cos0") {ff-^\n), (35) 

and observe that (^, 0) is a polynomial in even powers of ^ with the highest power given by 

^max^^-j _ 2 f m - 1 for m odd 



k=jn-i \m — 2 for m even. 

The conclusion in Eq. (36) completes our proof of polynomiality, Eq. (17), for {E'^ — E'^){x,^,t). 
Noteworthy, for an even moment the highest power could be m according to Eq. (17), but it is 
m — 2 in the model. In the next section we will see that this is a consequence of general large- A^c 
counting rules. 



6 The i:>-term 

The double distributions (DDs) of Refs. [1, 2] have the advantage of satisfying the polynomiality 
condition automatically. However, the description of non-forward nucleon matrix elements of two- 
body operators on the light cone in terms of DDs is incomplete, if one assumes them to be regular 
functions ^. This becomes manifest in the absence of the highest power in even (singlet) moments 
m of H'^{x,^,t) and E'^{x,£^,t). In other words, the coefficients hm™'^ and Cm™^ in Eq. (18) vanish. 
The description can, however, be completed by introducing the D-term D'^(z), where z = xj^ 
[21]. The D-term satisfies D'^iz) = —D'^{—z), has a non-vanishing support for \z\ < 1 only and is 
uniquely defined in terms of its moments [21, 22, 23] 

r dzz"'-^Di{z) = h'"^"'^ = -'''^"'^ for even m ^37^ 
J-i for odd m, 

^In principle it is possible to obtain a complete description (and to avoid the introduction of the D-term) by using 
an alternative definition of double distributions [22]. 
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with /im™^ and €■771 ^ defined in Eqs. (16, 17, 18). 

The D-term exhibits the following large- A^c behaviour [7] 

{D^ + D'^){z)=N^fiz), 

{D--D'){z)=N,f{z), (38) 

where the functions f{z) = 0{N^) are constants in the limit of large (and of course different for 
the different flavour combinations). It is convenient to expand the D-term in terms of Gegenbauer 
polynomials C^{z) for u = ^, the solutions of the leading order ERBL-evolution equations [33, 34], 
which govern its dependence on the renormalization scale /j,, 



D^{z) = {l-z') Y: dlcl'\z). (39) 

fc=l,3,5,. 



The constant d\ is of particular interest, as it determines the behaviour of the Z)-term (and of 
OFDFs) in the asymptotic limit of renormalization scale ji ^ oo [7] ^ . 

In the literature the following phenomenological approach for modeling became popular. One 
models the regular part of DDs (relying on educated guesses, physical intuition) and supplements 
the description by adding the Z)-term (using often information from the xQSM [24]). This approach 
yields DDs/OFDFs consistent with all general requirements including polynomiality. Many recent 
estimates and predictions for exclusive processes have been done in this way [24, 35, 36, 37]. 

The phenomenological interest motivates further studies in the model. Let us discuss first the 
fiavour-non-singlet combination (D" — D'^){z). From Eqs. (35, 36) we see 



dz z'^-\D'^ -D^){z)={) V m , (40) 

which means that also (D" — D'^){z) = 0. Note, however, that this is a zero at the order 0{Nf). 
We obtained this result from moments of {E'^ — E'^){x,^,t) = 0{N^) according to Eq. (10). Due 
to the relations (18, 37) it is (iJ" - H'^){x,^,t) = 0{Nc) which dictates the large-iVc behaviour of 
the flavour-non-sing let D-tcrm in Eq. (38), and not (^" - E'^){x,C,t) = O(N^). Consequently, the 
D-term does not contribute to (£"" — E'^){x,^,t) at the order of Nc considered here, which is just 
our observation in Eq. (40). The result (40) is thus a further demonstration of the consistency of 
the model, which respects the general large- A?c counting rules. 

From Eq. (27) we read off the model expressions for moments of the flavour-singlet {D^+D'^){z), 



1 



-1 



dz z'^-\D'' + D'^){z) = 2NcM^ (n| tV C>("*V) > 

n, occ 

•771 ^Tl— 1 

OM ^ 1_ Y {'^J^){fy\^rPmicose)ifr-^-\ meven. (41) 



j=0 

Recalling that Mn = 0{Nc), we see that all moments of {D'^ + D'^){z) are 0{N^) in the large-A^c 
limit, in agreement with the general counting rules in Eq. (38) also in this case. 

In principle Eq. (41) enables us to evaluate all moments of the flavour-singlet L>-term. Here 
we shall restrict ourselves to the evaluation of the second, i.e. lowest non-vanishing, moment. It is 



^D'^{z) mixes under evolution with the corresponding gluon D-term D^{z). Thus, in the flavour singlet case, it is 
strictly speaking ^ d\-\-d\, which is scale independent and determines the asymptotic behaviour of the Z)-term [7]. 
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related by J^^ dzzD^{z) = ^d\io the first coefficient d\ in the Gegenbauer expansion in Eq. (39). 
The expression for the flavour singlet d\ = d^ + df follows from Eq. (41) for m = 2 and reads 

di = -j N,M^ (^1 [f , |X| P2(cos e)\\n) . (42) 

n, occ ^ 

We evaluate di by means of the so-called interpolation formula [18], which consists in exactly evalu- 
ating the contribution to di in Eq. (42) from the discrete level (i!f^, and estimating the contribution 
from the negative Dirac continuum, dj""*, as follows. One rewrites d™'^* in terms of the Feynman 
quark propagator in the static background pion field J7, expands it in a series of gradients of the 
[/-field, and evaluates the leading contribution. The interpolation formula yields exact results in 
three limiting cases: (i) low momenta, |V[/| <C M, (ii) large momenta, |V?7| S> M, (iii) any mo- 
menta but small pion field, | log[/| ^ 1. One can thus expect that it yields reliable estimates also 
for the general case. Indeed, it has been observed that estimates obtained in this way agree with 
results from exact and much more involved numerical calculations within (10 — 20)% [18, 19, 20]. 

The expectation value of the operator in Eq. (42) vanishes in the grand-spin K = Q and positive 
parity state, which characterizes the discrete level, i.e. we have 

d^ = . (43) 

For the contribution of d™'^' we obtain 

where trp denotes the trace over flavour. In the calculation leading to Eq. (44) we took the (generally 
momentum dependent) mass M constant for simplicity. This simplification is legitimate for finite 
(or at most only logarithmically divergent) quantities, which do not (or only weakly do) depend on 
regularization. 

For constant M the continuum contribution df'^^ is logarithmically divergent and has to be 
regularized. The same type of divergence appears in the model expression for the pion decay 
constant /tt. Regularizing both quantities consistently in the same scheme allows to eliminate 
regulator and cutoff in favour of f-j^. It should be noted that further terms, containing three or 
more gradients of the [/-field and neglected in Eq. (44), are finite ^. 

In the hedgehog-ansatz the [/(x)-field is specified by the soliton profile P(r) with r = |x|. We 
take for P{r) the analytical arctan-profile 

P{r) = -2arctan-2- , R = M^^ = (350 MeV)"^ , (45) 

which is known to well approximate the self-consistent profile. In Eq. (45) R is the so-called size 
of the soliton. This enables us to continue the calculation analytically, and yields 

^cont = _ 3 7^2^ f^M^R^ . (46) 

With = QlMeV and the value Mn = 1170 MeV obtained in the xQSM with the analytical 
arctan-profile (45) - and "adding up" the zero of Eq. (43) - we obtain 

di = -9.46 at fj? ~ (600 MeV)^, (47) 
''For more technical details on the interpolation formula, and examples of similar calculations see Refs. [18, 19]. 
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the low scale of the model. In order to understand how this result can be used for phenomenology, 
we write down the entire sum rule for the second moment of (i?" + H'^){x,^,t) which reads 

1 

Jdxx {H^ + H") (x, 0) = Mq + ^ di e , (48) 

-1 

where Mq is the fraction of the nucleon momentum carried by quarks and antiquarks. In the xQSM 
at its low scale Mq is exactly unity since there are no gluons in the model and the entire nucleon 
momentum is carried by quarks and antiquarks [18]. Considering ratios of model quantities to be 
more reliable than absolute numbers, the prediction of the xQSM is 

-9.46 at ~ (600 MeV)^. (49) 



Mq 



Neglecting the scale dependence of the ratio di /Mq (and this we can do since the effects of evolution 
are smaller than the intrinsic accuracy of the model result in Eq. (49)) wc obtain with Mq ~ 0.5 
at a few GeV^ the prediction from the xQSM for experimentally relevant scales 

di ~ -4.7 at few GeV^. (50) 

The result in Eq. (50) agrees well with the number di ~ —4.0 obtained from the numerical results 
on {H^ + H'^){x,^,t) of Ref. [14] and quoted in Ref. [24] ^. The difference between the two numbers 
is due to the use of different regularizations. In Ref. [14] the momentum dependent mass M{p) 
was used as a regulator. Here we took constant M (and applied some regularization scheme, the 
dependence on which has been then removed in favour of /,r)- 

The advantage of our method is that we are in a position to exactly evaluate, at least in 
principle, any moment of the D-teim in the model. In Ref. [24, 38] this was not possible, not even 
in principle, due to limitations set by numerical accuracy. 



7 Summary and conclusions 

We have presented explicit proofs that the xQSM expressions for the leading flavour combinations 
in the large-iVc limit of the unpolarized OFDFs (i?" + i/'^)(x, ^, t) and {E'" - E'^){x, ^, t) satisfy the 
polynomiality condition. The method can straightforwardly be generalized to the case of flavour 
combinations sublcading in iVc, and to the case of helicity dependent OFDFs. The proof is an 
important contribution to the demonstration of the theoretical consistency of the description of 
OFDFs in the framework of the xQSM, and increases the credibility and reliability of theoretical 
predictions made on the basis of the xQSM. 

As a byproduct of the proof, we derived explicit model expressions for the moments of the 
flavour-singlet D-term, which are of phenomenological and theoretical interest. The explicit model 
expressions not only shed some light on the chiral dynamics underlying the D-term, but also 
enable exact model calculations of the coefficients in the Gegenbauer expansion of the D-term. 
We have demonstrated it by analytically computing the first coefficient di for the flavour-singlet 
I?-term reproducing the numerically estimated value previously obtained. Calculations of higher 
Gegenbauer coefficients are in progress. 

The Gs-symmetry transformation used in App. C is the "non-standard" time reversal transfor- 
mation [39] for quark fields in a static background pion field. We became aware of that after this 

^By numerically taking moments of + H'^){x,£^,t) from Ref. [14] and extrapolating to t = the value 
di/MQ ^ —8.0 has been found at the low scale of the model [38]. Neglecting the effects of evolution and using the 
phenomenological value Mq w 0.5 in Ref. [24] it was estimated di m —4.0 at experimentally relevant scales of few 
GeV^ in Ref. [24]. 
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work was completed thanks to Ref. [40]. Thus, as in QCD, the proof of the polynomiaUty property 
of OFDFs in the chiral quark-sohton model also relies on the principles of time reversal invariance 
(and hermiticity) and Lorentz-invariance, with the latter reduced to rotational invariance, once the 
soliton rest frame has been chosen to evaluate the model expressions for OFDFs. More precisely, 
due to the specific hedgehog ansatz for the soliton field, the model expressions are invariant under 
simultaneous rotations in space and flavour space. It is not surprizing that the deep connection be- 
tween flavour, space-time and chiral symmetry becomes apparent most easily here, in a relativistic 
chiral model of the nucleon based on the large- AT^, limit. 

Wc would like to thank M. V. Polyakov and M. Schwamb for discussions. This work has partly been performed 
under the contract HPRN-CT-2000-00130 of the European Commission. 

A Mellin moments of (iJ" + iJ^) and - E"^) 

H{x,^,t) = (ii"" + H'^){x,^,t) and E{x,^,t) = {E'' - E'^){x,^,t), Eqs. (12, 13), can be written 
compactly as 

A{x,C,t) = c^ E /^e-"(^^-^")/„(r^,zO,A) (51) 



with 



Ch = NcM^ , Ce = {2/3)iNcM^/{A-^f , 

Fff = (1 + 7°7^) , r^; = (1 + j^'j^Xt X Af , 



(52) 



and /n(r, A) defined as 



/„(^,z^A) = Jd 



'X e^^^$;(X + ^e^) r$„(X - ^e^) 



(n|r exp{—i—p^) exp(iAX) exp(—i—p^)\n) 



(53) 



where in the second line we made use of $„(X— a) = (X— a|n) = (X| exp(— ■jap)|n) and analogously 
for $;(X - a). For the m}^ moment M^\^,t) we obtain 



1 f jzo 



M^r\i,t) ^ Uxx^-^ A{x,i,t) = cj, /d^^""'/^ 



{xM-^-En) 



-'"n nocc 



-'"n nocc 



E 



d£^ 

27r 
d 



,m— 1 



d{iz^ 



Q\m-1 



dy 



e-^^°^"/n(rA,^°,A) 



m— 1 



g(_j^O)m-l 
m— 1 

E Er^-')^n-^- 



N n occ j^=0 



;0=0 

d{-iz^f 



zO=0 



(54) 



In the second line of Eq. (54), after the substitution y = xM^, the new integration limits [—M^, M^] 
have been replaced by [— cxd, oo]. This step is justified in the large-iVc limit where Mn = 0{Nc). 
The differentiation of /n(r, A) in the last line of Eq. (54) yields 



9^/„(r,z",A) 

d{-iz^)>' 



^0=0 j=0 



k 

^ E ( j ) ify exp(iAX) {f)''-^\n) . 



(55) 
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Inserting Eq. (55) into Eq. (54) we obtain 

Mr(C,0 = ^ E E("^fc"')%^E(J)HrA(p^)^exp(zAX)(p3)^-^|n), (56) 

N n,occ ^ = i=0 

which is the result quoted in Eq. (19) and in Eq. (29), respectively. 

B Analytical continuation to t = 

The moments Mj^(^,t) and Me(^, t) in Eq. (56) depend on the variables ^ and t only through 
the 3-momentum transfer A in the "large- A'^c kinematics" of Eq. (9). As only ^ and t are relevant 
variables, A can be expressed through ^, t and some arbitrary real angle a as 

A = (^-sina^-i- (2^Mn)2 , cosa^-f - (2^Mn)2 , - (2^Mn)^ . (57) 

The arbitrariness of a reflects the azimuthal symmetry of the physical situation around the 3-axis 
chosen for the space direction of the light cone in Eq. (9) . Since there is no real dependence on this 
angle, we can choose a particular value for a, or equivalently average over it. 

(i) (H"^ + H'^). To continue the moments Mi{{C,t) to t = 0, wc have to continue analytically the 
function FH{S,,t) = exp(iAX) in Eq. (19). Taking A as defined in Eq. (57) we obtain a function 
Fffi^jt) depending in addition on the angle a, and we remove the dependence on this arbitrary 
angle by averaging over a G [0,27r]. Taking X = |X|(sin0 cos(/>, sin 9 sm(f), cos 9) in spherical 
coordinates and expanding exp(iAX) in a series, we obtain for Fff(^,i) 

°° f-zlXlV 1 / / 

FH{^,t) = exp(iAX) = J2 II 2^1 V^~^ ~ (2C^n)2 sin^sin(a -((>) + {2^M^) cos 9 







(58) 

We see that i<//(^,t) is an analytic function of the variables ^ and t everywhere except for the 
twofold branching point t = — (2^Mi,j)^. Taking the limit t — > within a definite Riemann sheet, 
we obtain 

oo I '^^ I 

lim FH{^,t) = y (-^^^^nIXI) 1 a ^ ■ sin9sm{a - (^) +cosg') (59) 

alytical l\ ZTTJ \ J 



analytical 
itinua 
t^O 



continuation 1=0 Q 



with the plus and minus signs referring to the two Riemann sheets. The final result is real and 
does not depend on the choice of sheet 

_^ 27r ' 1 ^ ' 

— Jda sin 9 sm{a — (f>) + cos 9^ = — Jda sin 9 sin a + cos 9^ = Pi (cos 9) . (60) 



One arrives at the intermediate result in Eq. (60) by using the periodicity of the sine and cosine 
functions. The last equality is the integral representation of Laplace and Mehler for Legendre 
polynomials [41]. Inserting Eq. (60) in Eq. (59) we arrive at the result quoted in Eq. (20), 

lim ^ F^i^,t) = E ^"^'^^f^'""'^' Pl{cos9) . (61) 

analytical I \ 

continuation 1=0 
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(ii) - E"^). From Eq. (29) we see that the function Fe{^, t) - which is to be continued to the 
point i = in that case - reads 

_ iA^exp(iAX) _ , FH{i.t) 

= -t - (2^M,)2 (^2) 

with FH{i,t) from Eq. (58). In Eq. (62) V"^ = For ^ 7^ the hmit t yields 

hm , Fe{^, t) = V^f: (-i^^M \^\y~" p,(cos e) , (63) 

analytical ' ' I [ 

continuation 

where we took the hmit independently on numerator and denominator, since both limits exist, and 
used the result (61). The / = addendum in the sum in Eq. (63) is a constant and vanishes upon 
differentiation, so the sum starts with 1 = 1. In the next Appendix we will see that the addendum 
for Z = 1, which is singular at ^ = 0, does not contribute. Thus the moments M^\^,Q) are well 
defined for all ^. 

Using the identity ^ F{X.) = i[p,F(X)] in Eq. (63) we obtain the result quoted in Eq. (30). 

C Symmetries of the model 

Gg-symmetry. Consider the unitary matrix G5 with the following property 

GsT'^Gg 1 = (7'^)^ , Gsr^Gs 1 = -(r")^ . (64) 

In the standard representation of 7- and r-matrices G5 = 7-'^7^r^. Noting that V* = — (V)-^ holds 
formally, one finds that G5 transforms the Hamiltonian Hcs and the single quark states, Eqs. (2, 5), 
in coordinate-space representation as G^Hc^G^^ = H^^ and G5<I>„(x) = 'I>*(x) [17]. 

Let F be some matrix in Dirac- and flavour-space and define F = (G5FG5 ^)^. Let -F(X) be a 
function of X. Then 

{n\r{fyF{-X){p^r\n) = {-iy+'^{n\r{p^rF{-X){fy\n) . (65) 



Parity. The parity operator 11 = fl^^ = ^^P, where VF(X.)'P '^ = -F(— X), acts on the single 
quark states as Illn) = ±|n) according to the parity of the single quark state in Eq. (5). With the 
notation of Eq. (65) upon application of the parity transformation we have 

(n|F(^3)'F(X)(p3)"^|n) = (-l)'+'"(n| (7°F7°) {p^f (pF(X)p-i) {fr\n) . (66) 

(i) (H^' + H'^). Here F = (I + 7V) and f = (I-7V), and F(X) = \±f-Pi, where Pi = Picose) 
is set for brevity. According to Eq. (65) 

k k 

^ (J) (n|(i + ^W){fym'Pi{ff-'\n) = (-1)' E ( j) ("1(1 - 7V) [ff-'m'Pi {fy\n) . 

j=0 j=0 

(67) 

From Eq. (67) it follows that the contribution of the "1" in (1 -|- 7''7^) vanishes for odd k, while 
the contribution of "7*^73" vanishes for even k. Using (7*^7^)*^ = 1 for even k, and (7*^7^)^ = 'y^'y^ 
for odd k, this conclusion can be written compactly as 

k k 

X: (J) (n|(l + j'^'){fy\X\^Pi {f)'-^\n) = (J) (n|(7V)'(^')^|X|'P, {f)'-^\n) . (68) 
j=o j=0 
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Applying Eq. (66) to the intermediate result (68) yields 

k k 

E ( j) (^1 (tV)' ifVm'Pl {p^)'-'\n) = {-ly E (i ) (^1 (tV)' {py\^\'Pl {p')'-^n) , (69) 

j=0 3=0 

i.e. the expression vanishes if I is odd. The result (69) yields the expression in Eq. (21). 

(ii) (E" - £"^). Here T = (1 + 7°73)r" with f = (-1 + 7°73)r". It is convenient to take F(X) as 
in Eq. (63), F(X) = (^|X|'Pi) = (V^|X|'Pi). Note that the derivative acts on \±\^Pi only. The 
symmetry transformations (65, 66) yield 

k 

E (i) (n|(7V)'+'r»(p3)^(v''|X|'P,)(p3)'=->) 

3=0 

k 

= (-1)' E (^) (n|(7V)'+V« (p3)i(v''|X|'P,)(p3)*^-i|n) , (70) 

3=0 

i.e. also here contributions of odd / vanish. Applying Eq. (70) to Mg"^(^,0) and using the identity 
V^P(X) = ^[p^F(X)] we obtain the result quoted in Eq. (31). 

D Irreducible tensor operators 

The grand-spin operator is defined as 

K = L + S + t, (71) 

where L = r x p is the orbital angular momentum operator, S = ^757*^7 is the spin operator 
and T = ir is the isospin operator. K is the generator of simultaneous rotations R in space and 
SU (2)-flavour space. Let n be a unit vector defining an axis in space and flavour space and let a 
be an angle. The rotation R is then characterized by the unitary operator U {R) 

C/(i?) = ?7(n,a) =exp(-ian-K) . (72) 

The unitary transformation in Eq. (72) transforms single quark states into "rotated quark states" 
according to 

K 

U{R)\En,7r,K,M)= ^ \En,7r, K, M')D^^}j^iR) . 

M'=-K 

Here D^^)j^{R) denote finite-rotation Wigner matrices. 

The operators with integers L, M G [— L, L] are said to be irreducible tensor operators of 
rank L, if they transform into each other under the action of the unitary transformation in Eq. (72) 
according to 

L 

U{R)f^U\R)= E D^M'M{R)Th'- (73) 

M'=-L 

By considering infinitesimal rotations one arrives at a definition, which is equivalent to Eq. (73), 



L{L + 1) - M{M ± 1) f^±i , where = ± iK^ , (74) 
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but allows more easily to check whether a given operator is an irreducible tensor operator or not. 
The Hamiltonian H^g transforms as U(R)HcfiU^^{R) = HeS, since it commutes with K. This 
means that H^s is an irreducible tensor operator of rank zero. 

Further examples of rank zero operators are |X|, r • X or 7^7 • p, which satisfy Eqs. (74) for 
L = M = 0. The operators 



1 _ ^^(=FP^ -ij3^)/\/2 forM = ±l 



^^"^^ \if forM = 0, 

1 1 _ r i(Tr^ - iT^)/V2 for M = ±1 

l-'-blM 



ir^ for M = 0, 

|i ^ j i{fX^ -iX^)/V2 forM = ±l 
iX^ for M = 0, 

1 _ / ^(T7°7^ - i7°7^)/V2 for M = ±1 



d\M I .^0^3 fo, ^ ^ 0, 

rfii _ r i(T(T X p)i - z(t X p)2)/V2 forM = ±l . . 

^ " ^i(Txp)3 forM = 



are irreducible tensor operators of rank 1 since they satisfy Eqs. (74) for L = 1. 

Finally Legcndrc polynomials -Pi(cos 6) arc irreducible spherical tensor operators T/^ with m = 0, 
smce P^cos^) oc Y^{n)\ 

m=o = T'^|m=o- The spherical harmonics Y^(Cl) are irreducible tensor 
operators of rank I, since they satisfy Eqs. (74). 
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